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Spatial mapping of ices in the Oph-F core 
A direct measurement of CO depletion and the formation of C0 2 
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ABSTRACT 



Aims. Ices in dense star-forming cores contain the bulk of volatile molecules apart from H 2 and thus represent a large fraction of dark cloud 
I chemistry budget. Mm observations of gas provide indirect evidence for significant freeze-out of CO in the densest cores. To directly constrain 
the freeze-out profile of CO, the formation route of C0 2 and the carrier of the 6.8 yum band, the spatial distribution of the CO/C0 2 ice system 
\Q " and the 6.8 yum band carrier are measured in a nearby dense core. 

C^- • Methods. VLT-ISAAC, ISOCAM-CVF and Spitzer-IRS archival mid-infrared (3-20 ^m) spectroscopy of young stellar objects is used to 
■ construct a map of the abundances of CO and C0 2 ices in the Oph-F star-forming core, probing core radii from 2 x 10 3 to 14 x 10 3 AU or 
densities from 5 x 10 4 to 5 x 10 5 cm 3 with a resolution of ~ 3000 AU. 

Results. The line-of-sight averaged abundances relative to water ice of both CO and C0 2 ices increase monotonously with decreasing 
distance to the core center. The map traces the shape of the CO abundance profile between freeze-out ratios of 5-60% and shows that 
the C0 2 ice abundance increases by a factor of 2 as the CO freezes out. It is suggested that this indicates a formation route of C0 2 on a 
CO ice surface to produce a C0 2 component dilute in CO ice, in addition to a fraction of the C0 2 formed at lower densities along with 
Oh the water ice mantle. It is predicted that the C0 2 bending mode band profile should reflect a high CO:C0 2 number ratio in the densest 
q parts of dark clouds. In contrast to CO and C0 2 , the abundance of the carrier of the 6.8yimr band remains relatively constant throughout 
*_i the core. A simple freeze-out model of the CO abundance profile is used to estimate the binding energy of CO on a CO ice surface to 8 14 ±30 K. 
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^ , 1. Introduction IPontoppidan et alJ J2004J) introduced the technique of ice 

mapping at high (~ 1000 AU) spatial resolution - comparable 

Theory has long predicted that molecules freeze out onto dust t0 that of gas-phase maps. The procedure was demonstrated by 

grains in dense molecular clouds causing gas-phase abun- constructing a map of the distribution of water and methanol 

dances to drop by orders of magnitude. A series of recent mea- i ces in the outer envelope of the class protostar SMM 4 in 

surements of the gas-phase abundances of volatile molecules the Serpens star-forming cloud. The ice map of SMM 4 was 

such as CO and N?H + in dense cores have corroborated use d to show that the water ice abundace remained constant 

this conjecture (e.g . |Caselli et alj| 19991 |Bacmann et al||2002|; over a relatively wide range of densities, while the methanol 

iTafalla et alJl2004. CO ice has also been observed directly abundance increased sharply by at least a factor of ten within 

along a growing sample of isolated lines of sight toward both 10 000 AU of the center of SMM 4. In principle, ices are best 

embedded young stellar objects as well as background stars, mapped toward field stars located behind the cloud. However, 

showing that the CO mol ecules taken from the gas-phase re- such background stars are typically extremely faint in the mid- 

appea rs in the solid phase JChiar et all 19941 IPontoppidan et alJ infrared wavelength region. Therefore, disk sources embedded 



120031). CO ice is an excellent tracer of freeze-out processes be- m their parent molecular cloud core offer a convenient infrared 
cause it is the only abundant ice species known to form initially continuum against which the ices in the core can be mapped, 
in the gas-phase before adsorbing to a grain surface. w i tn appropriate caveats on the interpretation of the derived ice 

abundances. 
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Direct mapping of CO ice abundances is highly com- 
plementary to mapping of gas phase abundances, because 
ice maps are sensitive to depletion fractions as low as 5%, 
while gas-phase abundance maps trace high de pletion factors 
(n(ice)/«(gas)) of > 50% dCaselli et all fl 999). Furthermore, 
the depletion fraction as a function of gas density in a core 
is strongly dependent on the properties of the dust grain sur- 
faces. Direct maps of CO ice in a dense core therefore enables 
an independent measure of surface binding energies. 

This letter presents the first spatial map of CO and CO2 
ices in a dense molecular core. It will be demonstrated that by 
directly measuring the relation between ice abundances and 
gas densities, the following quantities can be constrained: I ) 
the dependence of CO depletion on gas density with a large 
dynamical range, 2) the CO-CO binding energy, 3) the for- 
mation route of CO2 and the 6.8 pm band carrier. The case 
stud y is the F core in the Ophiuchus molecular cloud complex 
(Mot te et alJ (11998)). Embedded within this region, are at least 
8-10 infrared-bright young stars. The lines of sight toward 5 of 
these young stars are used to probe the ices in front of each star 
in order to obtain a cross section of ice abundances across the 
minor axis of the core. 

2. Observations 

The individu al spectra are ta k en fro m Ivan Dishoeck et ail 
J2003I) and IPontoppid an et alJ d2003l) of the 3.08 pm and 
4.67yum stretching modes of solid H2O and CO, obtained 
with the Infrared Spectrometer And Array Camera (ISAAC) 
on the Very Large Telescope (VLT) The solid CO2 com- 
ponent is probed along the same lines of sight using the 
InfraRed Spectrograph (IRS) on the Spitze r Space T elescope 
(APR IDs 0009346048 an d 0009829888 2 ) JPontonnidan et al] 
l2005t lLahuis etaD 120061) as well as 5-16.3 pm spectra ob- 
tained with ISOCAM- CVF (TDTs 29601715 and 29601813) 
( [Alexander et alJl2003l) 3 . The Spitzer spe ctra were taken a s 
part of the Cores to Disks Legacy program fevans et all20 03). 
Combining these facilities, high quality 3-16jum spectra are 
available for 5 sources within a radius of 15 000AU from the 
center of the pre-stellar core, defined to be the 850 u m peak of 
Oph-F MM2 (a=\6 h 21 24.3, <5=-24°40 35, J2000) dMotte et alJ 
1998). Only IRS 46 has no suitable spectrum between 5 and 
lOjum. 

To determine the optical depths of the various ice bands 
continua were fitted using low-order polynomia. The CO and 
CO2 bands are sufficiently narrow to make the continuum de- 
termination relatively straight-forward and unbiased. For the 
3.08 pm water ice band a continuum is fitted to points between 

1 Partly based on observations obtained at the European Southern 
Observatory, Paranal, Chile, within the observing programs 164.1- 
0605 and69.C-0441. 

2 This work is based in part on archival data obtained with the 
Spitzer Space Telescope, which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology under a contract with 
NASA. 

3 Partly based on observations with ISO, an ESA project with in- 
struments funded by ESA Member States (especially the PI countries: 
France, Germany, the Netherlands and the United Kingdom.) 
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Fig. 2. The locations (triangles) of the sources used to construct 
the ice map plotted on the SCUBA 850 pm map. The star sym- 
bol indicates the center of the core. 

3.8 to 4.0 pm and a Zf-band photometric point from the 2MASS 
catalogue. Finally, a local continuum between 5.5 and E.Qpm is 
used to extract optical depth spectra of the 6.8 pm band. In con- 
sideration of the uncertainties in the continuum determination, 
care was taken to use the same "ice-free" regions and 2-order 
polynomia for all sources. Further discussion of the determi- 



nation of continua for extracting 


ice optical depth spectra can 


be found in e.g. 


Gerakines et al. 


dl995h: iDartois et alJ (Eo02): 


Keane et al. 


(2001 


. The ice spectra are shown on an 


optical 



depth scale in Fig. [2 while the locations of the sources relative 
to an 850 pm J CMT-SCUBA map obtained by the COMPLETE 
collaboration jRidee et all200 6j) are shown in Fig. [2 

The CO ice col u mn d ensities have been determined 
in IPontoppid an et alJ |2003). For the water 3.08 pm band 
and the CO2 15.2//m band, band strengths of 2.0 x 
10~ 16 cm mol ecule" 1 and 1.1 X 10~ 17 cm molecule 1 are used, 
respectively (Gerakines et alJll995l) . For the 6.8 pm band, a 
band strength of 4. 4 X 10~ 17 cm/ molecule" 1 is assumed, ap- 
propriate for NH| (Schu tte & Khannal f2003). The ice column 
densities are summarized in Tabled 

3. Abundance profiles of CO and C0 2 ice 

The basic ice map is constructed by determining the ice col- 
umn densities along each line of sight and ratioing with the 
water ice column density. The implicit assumption is that the 
3.08 jum water ice band traces the total H2 column density (see 
IPontoppidan et alJd2004l) for a discussion). This yields CO and 
CO2 line-of-sight averaged ice abundances relative to H2O ice 
as a function of projected distance to the center of the core. The 
radial map is shown in Fig. [3] 

It is seen that the abundances of pure CO, water-rich CO 
and CO2 all increase toward the center of the core, while the 
6.8 pm band carrier decreases slightly in abundance. The sharp 
rise in CO ice abundance is expected for CO freezing out from 
the gas-phase at high densities in the central parts of the core. 
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Fig. 1. Overview of the H2O, CO, 6.8 jjm and CO2 ice bands (ordered left to right) used to construct the ice map of the Oph-F 
core. The 3.08 fim water ice bands shown on the left-most panel are ordered according to optical depth, for clarity, as indicated 
in the panel. In the remaining panels, the sources are ordered according to projected distance to the core center, with the source 
furthest from the center (IRS 46) at the top. No 5-8 fim spectrum is available for IRS 46. Note also that the 6.8 fim band of IRS 
42 is filled in by emission likely due to PAHs. All the spectra have been shifted vertically, for clarity. 



Table 1. Line-of-sight averaged ice abundances in the Oph-F core 



CBRR 2422.8-3423 IRS 43 



IRS 42 



IRS 44 



IRS 46 



Distance to Oph-F MM2 [AU] 
N(H 2 0) [cm- 2 ] 
N(Pure CO)/N(H 2 0) 
N(CO in water)/N(H 2 0) 
N(C0 2 )/N(H 2 0) 
N(6.8)/N(H 2 0)> 



1.8e3 

3.9el8 

0.78 

0.19 

0.32 

0.08 

ISAAC/IRS 



4.3e3 

2.96el8 

0.32 

0.073 

0.24 

0.10 

ISAAC/ISOCAM 



8.2e3 

1.87el8 

0.15 

0.05 

0.184 

ISAAC/ISOCAM 



11.0e3 

3.04el8 

0.06 

0.04 

0.154 

0.12 

ISAAC/ISOCAM 



13.9e3 

1.17el8 

0.05 

0.0 

0.103 

ISAAC/IRS 



Assuming a constant water ice abundance of 9 x 1CT 5 relative to 
H2 in acc ordance with a model of th e CRBR 2422.8-3423 line 
of sight of Pontop pidan et ail ll2005l) . the average abundance of 
pure CO ice is seen to rise from 4.5 X 10~ 6 to 7.0 x 10~ 5 , and 
the total CO and CO2 abundances from 2.3 x 10" 5 to 12 X 10" 5 . 
This corresponds to a total CO depletion ranging from 12% to 
60%, assuming an initial CO gas-phase abundance of 2 x 10~ 4 . 

The increase by roughly a factor two in CO2 ice abundance 
is particularly interesting. Since CO2 molecules are formed on 
the grain surfaces, the increase in CO2 ice abundance along 
with the CO freeze out is a direct indication of a formation 
route associated with CO. A significant amount of CO2 is 
present in the outer parts of the Oph-F core where the densities 
are not yet high enough for the CO to freeze out. This can be 
interpreted as evidence for two distinct eras for CO2 ice forma- 
tion: The first takes place at roughly the same time as the bulk 
of the water ice is formed, which is known to happen as soon as 
the extinction, Ay, into the cloud reaches a specific, relatively 
low, threshold of roughly 3-5 magnitudes ftVhittetetalll 9881. 
This domain is not probed by the Oph-F ice map. The second 
takes place during the catastro phic freeze-out of CO th at occurs 
at densities of a few 10 cm ( J0rgensen et al . 2005) in which 
an almost pure CO ice mantle forms. This second phase of 
CO2 molecules should be easily detectable in the shape of the 



15.2 fim CO2 bending mode absorption band, which is highly 
sensitive to the molecular environment of the CO2 molecules. 
Specifically, the new CO2 molecules should be found to be di- 
lute in the CO ice with an average fraction of 1 :5, as determined 
by the relative CO and CO2 abundance increases in the Oph-F 
core. Since the CO freeze-out rate is a very strongly increas- 
ing function of density, it is expected that the C02:CO ratio 
will decrease in denser parts of the core. Thus, one would ex- 
pect to find a range of components in the CO2 bending mode 
in the centers of dense clouds. The spectral resolution of the 
ISOCAM-CVF spectra (A/AA ~ 50) is not sufficiently high to 
search for such band shape changes in the current data. 

Since the young stars that are used to estimate ice abun- 
dances in the core material are embedded in the core, signif- 
icant caveats apply. Heating of the core material may desorb 
CO ice within a radius of a few hundred AU of each source, de- 
pending on the luminosity. However, since the projected size of 
the core is 30 000 AU, desorption is likely to play only a minor 
role. A possible exception to this is IRS 46, which contains a 
signif icant component of warm molecular gas along the line of 
sight jLahuis et alJ2006l) . Also, the sources may be surrounded 
by remnant envelopes with different ice abundances than the 
surrounding core. The ice a bsorption toward CRBR 2422.8- 
3423 was shown bv lPontoppidan et iD§005) to be dominated 
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Distance to Oph-F MM2 [AU] 

Fig. 3. Radial map of CO and CO2 ices of the Oph-F core. Ice 
abundances relative to H2O ice toward young stars are plotted 
as functions of projected distance to the center of the core as de- 
scribed in the text. The CO ice has been split into a "pure" a nd 
"water-rich" component following lPontoppidan et all 12003). 

by cold core material, although a fraction of the water, CO2 
and 6.8 //m ices are likely located in the disk. However, it was 
found that the CO ice observed in this line of sight could not 
originate in the disk. 

4. Empirical determination of the CO-CO binding 
energy 

The CO depletion profile in Fig.[3]is modeled using a simplified 
freeze-out model, assuming a static core and that only thermal 
desorption occurs. For the purposes of this letter, the model pre- 
sented is intended as a demonstration of the method, rather than 
achieving very accurate results. The rate of ice mantle build-up 
is then given by: 

d«ice 

—\ — = ^ads - ^des, (1) 

at 

Where R^ s = Vq exp[-dH/kT] x «co,ice x /? is the desorp- 
tion rate and R Rds = «co, g as x "dust x nd 1 x -\J3kT/mco X / 
is the adsorption rate. f3 is a factor taking into account that 
CO only desorbs from the top monolayer (i.e. 1st order des- 
orption for sub-monolayer coverage and 0th order desorption 
for multilayers), d = 0.05 /urn is the grain radius, / = 1 is 
the sticking coefficient, vo is the frequency of the CO stretch- 
ing mode and «co,gas and «dust are the number densities of CO 
gas and dust particles. Solving for n; ce yields a time depen- 
dent ice density given a gas density and (assumed identical) 
gas and dust temperature. This rate equation reaches an equilib- 
rium on relatively short time scales (<< 10 6 years). Therefore, 
given a temperature and density structure, one can in principle 
solve for the binding energy of CO, dH. For the density pro- 
file of the core, a Bonnor-Ebert sphere with a central density of 
5. 5 x 10 5 cm" 3 is used with a temperature of 15 K as suggested 
bv lMotte et al.ldl998l) and references therein. This simple pro- 
cedure yields a CO on CO binding energy of 814 ± 30 K. The 
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10 5 10" 
Maximum gas density along LOS [cm -3 ] 

Fig. 4. Observed and modeled CO ice abundance for Oph-F. 
The curve shows the simple freeze-out model for the best- 
fitting CO binding energy of dH/k = 814K. The "solid CO" 
abundance is the sum of the pure CO, CO in water and the CO2 
ice excess abundances. This assumes that one CO molecule is 
used to produce each CO2 molecule. The error bars indicate a 
20% uncertainty. 

uncertainty reflects different results obtained if varying the cen- 
tral density and water ice abundance by 50%. Additionally, the 
derived binding energy scales linearly with the assumed dust 
temperature. Considering the simplifie d analysis, this result is 
consistent with that recently found by Bissch op et alJ J2OO6) 
from laboratory experiments. 

The radial map constitutes the first direct observation of 
the freeze-out profile of CO on dust grains in a prestellar core 
previously inferred indirectly from observations of millimetre 
lines of molecules. Additionally the observed increase in CO2 
ice abundance toward the center is the first quantitative obser- 
vational evidence of the formation of CO2 from CO on the sur- 
faces of dust grains. These are observations that would not have 
been possible with single lines of sight. 
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